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INTRODUCTION 


Orbiter operations involving Spartan and other free-flying objects need the ability to see and detect such 
objects optically. This ability depends primarily on the brightness of the object rclati’i ’Q to other sources pre- 
sent. The present analysis develops a general mathematical model by which realistic estimates of brightness 
and visible range can be made, based on a simplified model of the free-flyer as an optical reflector, its attitude 
progran, and the predicted ephemerides of the Orbiter and the sun. 


Because of the uncertainties in predicting actual Spartan-Orbiter separation distances, the visibility 
prediction problem is formulated here to compute maximum working range from the other variables. Work- 
ing range is defined as that separation distance within which visual or optical tracking is operationally feasible, 
e.g., within which the brightness is sufficient for reliable detection by the Orbiter star trackers (third stellar 
magnitude). 


The Spartan reflector model is illuminated both by direct solar radiation and radiation reflected by the 
earth (albedo), producing a model source of defined directional intensity. Its intensity in the Orbiter direction 
(along orbit) yields the desired maximum range directly. The required geometric and photometric calculations 
involve a number of angles in space, which are readily computed from the basic directions defining their sides. 
The time-dependent directions are determined by straightforward calculation from fundamental relationships 
and constants. 


In practical application, a suitable series of computed range predictions, covering all or part of a Spartan 
free flight, can be plotted against time and compared to a set of actual range predictions or envelopes. The 
results permit determination of the intervals of usefulness of optical tracking for the mission, either as a 
primary or auxiliary mode. They may also suggest consideration of remedial measures, such as improving 
Spartan reflective properties or limiting the separation distance of Spartan from Orbiter. 
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azimuth of coordinate plane. A , 3 is azimuth at axis 1 of plane of axes 1 and 3. 
range-intensity correlation constant 
stellar magnitude definition constant 
non-negative cosine 

fraction of day from O*’ UT (universal time) 
day number (of year) 

distance (range), Orbiter to Spartan, meters 
illumination, lux (lumens/m^) 
earth center to earth albedo element unit vector 



eclipse function: eel $s(o) ~ 


when cos $g( © ) 



free-flight elapsed time (Spartan) 

Spartan face normal (unit vector) 

Spartan altitude, meters 
orbit inclination, ® 

luminous intensity, candelas (lumens/steradian) 

Julian date (day number counting from Jan. 1, 4713 B.C.) 

Julian date of last Jan. 0 
longitude of sun 

secular component of longitude of sun 
mean anomaly of sun 
mission elapsed time 

orbital revolution number (Spartan released N = 1) 
earth albedo element to Spartan unit vector 
distance, earth albedo element to Spartan, meters 
radius of earth, 6378.14 km (equatorial) 
area, m^ 

surface area of 1 AU radius sphere, m^ 

Spartan position unit vector 

time in Julian centuries from the new standard epoch J2(XX).0 = 2000 Jan. 1.5 = JD245 1545.0 
(Ref. 4) 

time since last ascending node passage, sec 
Spartan negative velocity unit vector 
right ascension angle on celestial sphere, ° 
off-nadir angle of element E (from Spartan), ° 
declination angle on celestial sphere, ® 
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€ obliquity of ecliptic (angle to equator), ° 
f sunlight reflection angle from earth albedo element, ° 

7] sunlight incidence angle on Spartan face, ° 

6 Spartan angular displacement from ascending node 

X orbital longitude, measured from ascending node of orbit, ° 

sunlight reflection angle from Spartan face, ° 
p reflectance = 0.30; ^j’s given Figure 1) 
cr albedo flux incidence angle on Spartan face, ° 

T orbit period, sec 

total luminous flux of sun, 3.79 x lO^* lumens 
^ sunlight incidence angle on earth albedo element, ° 

^s(o) sun-earth-Spartan angle (“solar colatitude” of Spartan) 

\p geocentric angle between albedo integration element and Spartan radius 
fl right ascension of ascending node of orbit, ° 

<j3 azimuth of albedo integration element, measured CW from Spartan meridian 
A ascending node 

^ sun (from earth) unit vector 

Subscripts 
A albedo 

a dummy angle symbol 

D direct (solar illumination) 

d deployment (Spartan) 

E end of mission 

ES albedo integration element to Spartan line 

e earth center to earth albedo element line 
f Spartan face normal 

i index identifying Spartan faces (i = 1, 2, 3, . . .6) 

j index identifying circular boundaries in earth albedo zone integration (j = 0, 1 , 2, . . .7) 

k index identifying sectoral boundaries in earth albedo zone integration (k = 0, 1, 2, . . .17) 

n rotation serial number 

R recovery 

S Spartan 

V along velocity vector 

w dummy axis index 

1,2,3 floating axis identifiers 
3m third stellar magnitude 

o sun 

® earth 

3 


Superscripts 

• derivative wrt time 
unit vector 

w complement of (angle), a = 90° - a 

~ supplement of (angle), a = 180° -a 

conjugate (explement) of (angle), a = 360° - a 
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CONTINUOUS PREDICTION OF SPARTAN VISIBILITY FROM ORBITER 
OVER MODELED FREE-FLIGHT MISSION 


SUMMARY 

For a Spartan free-flight mission specified by: 

1) Spartan attitude program, 

2) Spartan optical reflective properties, 

3) Orbit elements, 

4) Time phasing of attitude, orbital, and solar positions, 

compute and plot vs. time the estimated range at which the apparent brightness of Spartan is equal to some 
minimum useful level (e.g., third stellar magnitude for Orbiter star trackers). Considering only direct solar il- 
lumination of Spartan, the third magnitude range dj^j^ is derived from the luminous intensity Ip of Spartan 
by 

‘^3mD ~ U (1) 

where C = 2484.37 m/lumen'/^ (see CONSTANTS) 

_ P. Sj cnn Tji cnn /ij eel 

7T 

1 = 1 


Including the additional illumination of Spartan by sunlight reflected from the earth (albedo), the third 
magnitude range becomes 


dam = Cflo + Ia)'/^ = C I'/^ 


( 3 ) 


where 


2Eo(e)/)^ 6 


EE enn^j,^ Sj cosfj cnnCTjjkPiSiCnn/tj rj 
j =o k = o 


( 4 ) 


The above symbols are identified under NOMENCLATURE. Derivations and values are given in the discus- 
sion which follows. 


CONSTANTS 

Range-Intensity Correlation Constant 

Equations (1) and (3) relate range (distance d 3 „) for third magnitude brightness (illumination E 3 „,) to the 
luminous intensity of the “source” (Spartan). They are derived from the stellar magnitude definition 

100.4m E = C„ 

and the geometric relationship between intensity and illumination 

E = I/d2 (6) 

C„, is a constant (arbitrary but established) which can be evaluated from a known reference, such as the sun’s 
luminous flux = 3.79 x 10^* lumens and designated magnitude m= —26.8 at earth. Substituting (5) in (6) 

d3^ = (100-4-"/C„)'/^ I‘/^ = C 1'/^ 
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where C = 100 - 2 (m-mQ) = 2484.37 m/lumen 

S^y is the surface area of the sun-centered spherical surface through the earth over which the sun’s total 
luminous flux is considered distributed in computing the flux density (illumination at the earth (1 AU 

= 1.496 X 10'' m. away). 

Spartan Reflector Model 

For the purposes of this analysis, the optical reflective properties of Spartan 1 are approximated by the 
Spartan Reflector Model shown in Figure 1 . The figure lists the approximate face areas Sj and diffuse reflec- 
tance values Pj for use in Equations (2) and (4). 

Earth Albedo 

The average earth albedo value = 0.30 is adopted for use in this analysis. This value, taken from 
Reference 1, appears to be conservative. The illumination of Spartan by earth albedo flux is estimated via the 
approximate integration scheme shown in Figure 2, which defines the area element Sj. 


ANGLES 

The six generic angles which appear in the two intensity equations ((2) and (4)) are angles between 1 1 
directions in space: the Spartan position vector"?, the negative velocity vector (-^), the sun (from earth center) 
o , the earth albedo element (from earth center)^, the albedo element-to-Spartan line R£g,and the six Spartan 
face normals fj. These directions are all derived in the next section and expressed in equatorial coordinates 
(right ascension a and declination 5 on the celestial sphere). The required angles between the individual pairs 
are then determined using a single general relationship (spherical law of cosines): 

v_/ V-/ w 

cos a = cos5j cos 52 + sin5j sin52 cos (02 * 0 ) 


Table 1 lists the six generic ray angles, specifies their defining directions (sides), and tabulates their coor- 
dinates for entry into Equation (7) to compute the size of the particular angle defined. 


Table 1. Ray Angle Determination 


Angle 

Initial Side (1) 

Terminal Side (2) 


Description 

Direction 

Coordinates 

Description 

Direction 

Coordinates 

a 



5, 


“2 

52 


/A. 

Spartan face normal, f- 



Sun line, © 


5o 


Spartan face normal, f|^ 

"fi 

5fi 

Negative velocity vector, — v 


5-v 


Albedo element normal, e" 



Sun line, 0 ^ 

“0 

5o 

1 

Spartan position vector, 


“s 

Sun line, u 




Albedo element normal, 'e" 

“e 


Element-to-Spartan line, 

^ES 

^ES 

f 

Spartan face normal, f. 



Element-to-Spartan line, 

^ES 

•^ES 


cos a 



^ ^ W W w , 

cos8, cosSj + sin5j sindj cos (qj - a,) 


(7) 
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Reflectance 


Face 

Normal 

Area, S, 

(diffuse), Pi 

1 

X 

1.301 m2 

0.05 

2 

-X 

1.301 m2 

0.1 

32 

3 

Y 

2.694 m2 



4 

-Y 

2.694 m2 



5 

Z 

2.357 m2 



6 

-Z 

2.357 m2 




Z 
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DIRECTIONS 


Sun (0) 

The coordinates of the sun at any instant can be computed using formulas from Reference 2: 


= arctan (cos e tan L) 

(8) 

= arcsin (sin e sin L) 

(9) 


where e = 23.441 ° 

L = X + l.°916sin M + 0. “^20 sin 2M 
where X = 280. °460 + 36000. °772 T 
M = 357.°528 + 35999.050 T 
where T = (JD - 2451545.0)736525 

The Julian date JD of a given instant can be determined in various ways. A convenient way to accommodate a 
range of launch dates is to express JD as the sum: 


Julian Date of preceding January 0 

JDny 


Day number (of year) 

DN 

= 

Fraction of day from O'’ UT 

DF 

= 

Julian Date (sum of above three) 

JD 



Spartan Position Vector (s) 

V-/ 

Referring to Figure 3, the Spartan equatorial coordinates and Sj are found by solving the right 
spherical triangle ^GL to obtain 


tan Xj = cos i tan $ 

(10) 

N-/ 

COS 65 = sin i sin 6 

(11) 


The desired Spartan coordinate is obtained by adding the right ascension of the orbit ascending node fl to 
the Spartan orbital longitude Xj 


«s = fio + [(N-1) T + t] Q + arctan (cos i tan 6) (12) 

where = right ascension of ascending node immediately preceding Spartan release 
N = orbital revolution number (Spartan released during revolution N = 1) 

T = orbit period 

t = time since last ascending node passage 
Q = nodal precession rate 
d = k = 360Vr 
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Figure 3. Geometric Relationships Determining Directions 


Negative Velocity Vector Direction (-v) 

The Spartan velocity vector direction at a given moment coincides with that of the radius vector 90° 
ahead of Spartan. Substituting Ov = 0 + 90° into Equations (10) and (11) 


tan Xy = -(cos i)/tan $ (13) 

cosSy = sin i cos 6 (14) 

Again the orbital longitude \ is converted to right ascension by adding Q, the right ascension of the node 

= Qq + [(N - 1)t + t] Q + arctan (cos i)/tan d (15) 

>A. 

Albedo Element Direction (e) 

Referring to Figure 3, the coordinates and 5^ of the albedo element E are found by solving the oblique 
spherical triangle GNE yielding 


cos = cosSj cos \p + sinSj sin\p cosco 
sin (cKj-as) = (sin^ sinoj)/ sin6g 


(16) 

(17) 


xp and o) are explicitly defined and used to specify the grid of elements for the summation approximating the 
total albedo illumination of Spartan (see Figure 2). 

Element-to-Spartan Direction (Res) 

The element-to-Spartan linejdirection Rgs coordinates can be obtained from Figure 3 by producing a line 
from the earth center parallel to Rgg. If the radius OE is swung 90° clockwise along the great circle of EG to 
T, it coincides with the radius OT, which is parallel to the surface tangent at E. If the radius is swung back 
through the elevation angle f to R, it coincides with the desired parallel to Rgg through O. The desired coor- 
dinates are those of R, ags and 655, found by solving the oblique spherical triangle RNG 


cos Sgs = cos8^ cos(90°-^-f) -1- sinSg sin(90° - ^-f) cos S 
sin (as - a^g) = sincocos(^ -I- f )/ sinS^g 

The zenith angle f , as well as the distance rj, can be found by solving the oblique plane triangle EOS 

fj = K + + h)2-2r^ (r^ h) cos 

sinf = (r^ + h) (sin^)/rj 

Spartan Face Normals (fj) 


(18) 

(19) 

( 20 ) 
( 21 ) 


The time-varying orientation in space of the Spartan face normals is determined by the Spartan attitude 
program, which is preplanned based on observing and operating requirements. The attitude program is 
typically specified by an initial orientation (direction of one axis and rotational position of Spartan about it), 
followed by a timed sequence of rotations about selected axes. The specification format is illustrated by Table 
2, which is the basic attitude program planned for Spartan 1 . 
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While a plan like that of Table 2 completely determines the Spartan attitude program, it does not give ex- 
plicitly the axis (face normal) directions needed here for photometric calculations. These can be computed us- 
ing Equations (22) through (27), following the cyclic substitution routine illustrated in Table 3. 


Table 2. Attitude Program for Spartan 1 (4-23-84) 


Time 


Rotation 


Orientation After Rotation 


Julian 
Date, JD 


Rev. No., 
N 


Time 
since 
Node, t 


Type 

No. 


Serial No. 
V 


Axis* 


Displace- 
ment, A A 


Pointing* Axis 


Rt. Asc. 
a 


Deci. 

a 


up* Axis 


Azim.** 

A 


(days) 


(sec) 


(deg) 


(deg) 


(deg) 


(deg) 


(To be 
deter- 
mined. 
See Fig- 
ure 6) 


1997 
2085 
2217 
3458 
3541 
4977 
5045 

5117 
X 5329 

55 
136 
2181 
2266 
3472 
3552 
4976 
5045 

5118 


2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 
17 


0 

1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 
17 


Roll 

Ptich 

Roll 

Pitch 

Roll 

Pitch 

Roll 

Pitch 

Roll 

Pitch 

Roll 

Pitch 

Roll 

Pitch 

Roll 

Pitch 

Roll 


0.0 

-28.547 

-175.866 

15.889 

-88.899 

15.388 

-59.593 

-23.464 

-38.016 

32.478 

-67.477 

0.276 

-108.397 

12.347 

-86.899 

15.091 

-61.523 

-22.860 


20.662 I 45.101 


265.559-*- Initial 
Orient. 


J ^ * 

three orientation coord- 
inates for all six axis 
directions ( + X, +Y, +Z, 
the face normals) to be 
computed using algorithm 
shown schematically in 
Table 3. 


•a 3 

cd o 

II 

'5 Xi 

u cS 

< CJ 


4, 5,. 

5, 6,. 


.13 

.14 


e 

o 

o 

Q 


5329 

55 

136 

2181 

2266 

3473 

3561 

4989 

5052 

5118 


18 

19 

20 
21 
22 

23 

24 

25 

26 
27 


18, 28,. 

19, 29,. 

20, 30,. 

21, 31.. 

22, 32.. 

23, 33,. 

24, 34,. 

25, 35,. 

26, 36,. 

27, 37.. 


.108 

.109 

.110 

.111 

.112 

.113 

.114 

.115 

.116 

.117 


Pitch 

Roll 

Pitch 

Roll 

Pitch 

Roll 

Pitch 

Roll 

Pitch 

Roll 


-38.016 

32.478 

-67.477 

0.276 

-108.399 

15.804 

-99.619 

15.284 

-48.344 

-22.852 




t n 


15, 18. 21 

16, 19, 22 


5327 

55 

136 


28 

29 

30 


118, 124, 130 

119, 125, 131 

120, 126, 132 


Pitch 

Roll 

Pitch 


-38.016 

32.478 

-67.477 


17, 20 


^ c 
c/3 o D, 


2214 

2294 

4151 


31 

32 

33 


121, 127 

122, 128 
123, 129 


Pitch 

Roll 

Pitch 


67.477 

-32.478 

38.016 






23 


2188 

2266 

2303 


34 

35 

36 


133 

134 

135 


Roll 

Pitch 

Roll 


14.662 

-22.366 

50.000 


< S U 


-Ticks are asc. node crossings 


*Roll and Pitch rotations are about the X (Pointing) and Y axes respectively, positive clockwise when looking toward the 
origin along the positive axis. See Figure 1. 

** Azimuth angles are measured in the YZ plane from north to the Z axis, positive clockwise when looking toward the origin 
along the positive X axis. 
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Table 3. Completion of Body Axes Coordinates Matrix via Cyclic Algorithm 
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present azimuth A„. 











cos 52 = sin Ajj cos5j 

( 22 ) 

sin5j 

cos A 21 = w 
sin 62 

(23) 

cos X 2 = -sin A , 3 cos A 2 , 

(24) 

V-/ V-/ 

cos 63 = cos A, 3 sin 6 , 

(25) 

\m/ 

cos 6 , 

sin A 32 = 

sin53 

(26) 

cos X 3 = cos Aj 3 sin A 32 

(27) 


The geometric basis of Equations (22) through (27) is illustrated in Figure 4. The main sketch depicts the 
overall geometric relationships and defines the pertinent variables. The two small sketches are details of 
quadrantal spherical triangles which appear on the main sketch. Their solutions via their polar triangles yield 
Equations (22) through (27). 

The additional relationships needed for formal completion of Table 3 are straightforward: 


q ;2 = <^1 + ^2 

ctj = o![ + Xj (29) 


a_^ = ± 180‘ 


5-w = -K 


w = X, Y, Z 


(30) 

(31) 


BODY REFERENCE AXES 

The body reference axis system used in this analysis (Figure 1) is chosen for convenience, clarity, and 
compatibility with conventional usage. It is referred to as the “Pointing” system because the positive X axis is 
in the experiment pointing direction. The Z axis is conventionally “up” from the REM base in the Orbiter. 
The Y axis is to the left or port (along X-locator rods) when facing in the pointing direction, yielding a right- 
handed system. Rotations are considered positive which appear clockwise when viewed toward the origin. 

The Pointing system is shown in Figure 5 in relation to physical features of Spartan, and also in relation 
to the “ACS” system and two standard Orbiter reference systems. Note that all directions in the ACS system 
are opposite to those of the Pointing System. Positive rotations have the same physical direction in both 
systems because opposite look directions in defining clockwise rotation offset the opposite axis directions. 

COMPUTING AND PLOTTING SCHEME 

Generating the required range vs. time plots is a stepwise process requiring an arbitrary number of 
discrete range calculations, depending on the length of the mission, the particular Spartan attitude program, 
etc. Selecting a suitable calculating and plotting procedure may require iteration in some cases, but for a 
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practical attitude program in the form of Table 2, the simple procedure illustrated in Figure 6 will be a convenient 
tentative choice. 

Figure 6 also illustrates the time phasing relationships between the Spartan attitude program, which is 
phased primarily to the nodal passages . . .), STS mission events, and the basic Julian Date time 

reference system. As noted in Table 2, the Julian Date phasing is conveniently left open in preflight modeling, 
primarily to accommodate likely variations in launch date and time. 


Reference Axis Systems 


Face 

Pointing 

ACS* 

Body AXES** 

Orbiter 

Structural 

Identifying Features 
Pointing (front) 

X 

(-X) 

X 

Body** 

-X 

Antipointing (rear) 

-X 

Roll (X) 

-X 

X 

X-Locator Rods (port) 

Y 

(-Y) 

-Y 

-Y 

SOB Side (starboard) 

-Y 

Pitch (Y) 

Y 

Y 

Grapple Fixture (up) 

Z 

(-Z) 

-Z 

Z 

REM (down) 

-Z 

Yaw(Z) 

Z 

-Z 


*The ACS axes are left-handed if taken in the roll, pitch, yaw order. The others are right-handed. 
“See Reference 3. 



Figure 5. Spartan Body Reference Axes 
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Time, 



Figure 6. STS/Spartan Events Time Line with Trial Visibility Range Plot Format 
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